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ABSTRACT 

We measure the clustering of non-quasar X-ray AGN at z = 0.7— 1.4 in the AEGIS field. Using the 
cross-correlation of 113 C/ionrfra-selected AGN, with a median log Lx = 42.8 erg s~^, with ~5,000 
DEEP2 galaxies, we find that the X-ray AGN are fit by a power law with a clustering scale length of 
ro = 5.95±0.90 Mpc and slope 7 = 1.66±0.22. X-ray AGN have a similar clustering amplitude as 
red, quiescent and 'green' transition galaxies at z ~ 1 and are significantly more clustered than blue, 
star-forming galaxies. The X-ray AGN clustering strength is primarily determined by the host galaxy 
color; AGN in red host galaxies are significantly more clustered than AGN in blue host galaxies, with 
a relative bias that is similar to that of red to blue DEEP2 galaxies. We detect no dependence of 
clustering on optical brightness. X-ray luminosity, or hardness ratio within the ranges probed here. 
We find evidence for galaxies hosting X-ray AGN to be more clustered than a sample of galaxies with 
matching joint optical color and magnitude distributions. This implies that galaxies hosting X-ray 
AGN are more likely to reside in groups and more massive dark matter halos than galaxies of the 
same color and luminosity without an X-ray AGN. In comparison to optically-selected quasars in the 
DEEP2 fields, we find that X-ray AGN at z ~ 1 are more clustered than optically-selected quasars 
(with a 2.6(7 significance) and therefore may reside in more massive dark matter halos. Our results 
are consistent with galaxies undergoing a quasar phase while in the blue cloud before settling on the 
red sequence with a lower- luminosity X-ray AGN, i/they are similar objects at different evolutionary 
stages. 

Subject headings: cosmology: large-scale structure of the universe — galaxies: active — galaxies: 
high-redshift — X-rays: galaxies 



1. INTRODUCTION 

It has recently become clear that black holes play a 
key role in the evolution of galaxies. It is now under- 
stood that most galaxies have su permassive black holes 
in their nuclei (for revie ws, see e.g. iRichstone et al.lll99a 
iFerrarese fc Fordll2005[ ) and that these black holes power 
active galactic nuclei (AGN). There are strong observed 
correlations between black hole mass and galaxy prop- 
erties such as the stell a r velocity dispersion in the bulge 
(jCebhardt et al.l 120001 : [Ferrarese fc MerrittI l2000( ). that 
indicate some form of feedback or connection between 
the growth of black holes and their parent galaxies. 

It remains unclear how galaxies and AGN co-evolve 
and what impact the AGN have on the evolution of their 
host galaxies. It is also not known what the accretion 
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mechanism is for AGN and what their fueling source is. 
Measuring the environments and clustering properties of 
AGN is key to understanding their nature: constraining 
how they are triggered and fueled, inferring their life- 
times, and determining which galaxy populations host 
different kinds of AGN (quasars, Seyferts, radio-loud 
AGN, etc.). Additionally, clustering measurements con- 
strain the masses of the dark matter halos which host 
AGN and allow various types of AGN to be placed in a 
cosmological context. These pieces are all necessary to 
understand the relevance of AGN for galaxy evolution. 

Cosmological numerical simulations, both hydrody- 
namic and semi-analytic, predict the co-evolution of 
black holes and galaxies as well as the lifetimes and 
fueling mechanisms for AGN. In these simulations 
the AGN physics itself is not well resolved and of- 
ten assumptions are made regarding the details of 
accretion and feedback from the AGN. Most mod- 
els of high luminosity AGN assume that mergers 
drive the AGN accretion (e.g., |K auffmann & HaehnelJ 
I2OOOI : ISpringel. Pi Matteo. fc Hernquistll2005af l. though 
lower-luminosity AGN may be fueled in other ways 
(|Hopkins fc Hernguistl [2006) . Recent theoretical mod- 
els often include some form of AGN feedback, moti- 
vated by the observational links between the properties 
of galaxies and their central black holes; this feedback 
can cause a decline in the star formation activity of 
massive galaxies, creating red ellipticals. Many mod- 
els find that AGN feedback of some sort is in fact re- 
quired to cr eate the red seque nce of quiescent galaxies 
(though see iNaab et aP ()2Cl07| ) and iDekel fc BirnboimI 
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(|2008f ) for counter examples). For massive galaxies, 
accretion of gas from cooling flows in dense environ- 
ments may produce relatively low-luminosity AGN that 
in turn heat the bulk of the cooling gas and prevent it 
from fallin g into the galaxy's center to f orm stars (e.g., 
[Cattane o et al.ll2006l : fCroton et al.ll2006( l. Other models 
propose that mergers trigger quasars that drive outflows 
that sweep away the gas, not allowing it to f orm stars 
(e.g. JSpringel. Pi Matteo. fc HernqlIistil2005aD . 

These theoretical models can be tested using observa- 
tions of the location of AGN within the web of large-scale 
structure. Measuring the clustering of both low and high 
luminosity AGN is crucial for determining if their fueling 
mechanisms differ and searching for links between AGN 
fueling and local environment. To measure clustering 
properties accurately requires large samples to overcome 
both Poisson errors and cosmic variance. 

It is now known that optical surveys fail to detect 
a large fraction of the full AGN popul ation, particu- 
larly w eak or obscured AGN (e.g., Mushotzkv ct al.! 
2000t [Alexander et all 120011: Pushotzky 2004; 



Brandt fc Hasingeil 12005 ). X-ray detection is generally 



recognized as a more robust way to obtain a relatively 
unbiased A GN sample, though it does not detect a ll 
AGN (e.g., iDonlev et all [20051 : 'Hcckman ct al.' ' 2001 
espec ially optically-varying AGN (Morokum a et all 
120081 ) . As soft X-ray sample s can miss obscured AGN 
(e.g.. IComastri fc Fiorel[20M ). an ideal sample includes 
objects detected in either the hard or soft bands; 
however, even the deepest of the current X-ray surveys 
are likely to miss the most heavily obscured " Compton 
thick" AGN (jGilh. Comastri. fc Hasin"geill2007[) . 

Most previous attempts to measure the clustering 
of X-ray AGN in ROSAT, XMM-Newton and Chan- 
dra surveys have been fundamentally limited by a lack 
of spectroscopic redshifts, resulting in conflicting val- 
ues for the inferred correlation amplitude. Most stud- 
ies measure the projected angular clustering of X- 
ray sources (lAkylas, Georgantopoulos. fc Plionii 120001: 
Yang et al.' '200 3: iBasilak os et al 120041 : iMivIn et al.l 
2007; Plionis ct al. 2008); to interpret the angular clus- 
tering the redshift distribution of the sources, dN/dz, 
must be known to fairly high precision, as the inferred 
correlation length is particularly sensitive to the range 
of redshifts spanned. Lacking a large sample of spectro- 
scopic redshifts for these sources, most authors assume 
a redshift distribution using the observed X-ray lumi- 
nosity function and AGN population synthesis models. 
However, errors in the dN/dz distribution are gener- 
ally not propagated to the inferred correlation length, 
so that the quoted errors on the clustering amplitude are 
highly underestimated and often the inferred correlation 
le ngths are extremely la rge (e.g., tq = 9 — 19 h"^ Mpc 
in lBasilakos et al.l (|2005f l). Additionally, most studies do 
not include cosmic variance errors, which can be large for 
the relatively small fields with deep X-ray data. 

Recently, several studies have employed spectroscopic 
redshifts t o measure the direc t correlation function of X- 
ray AGN. iMullis et all (|2004D use data from the ROSAT 
North Ecliptic Pole survey to measure the clustering of 
219 bright {L^ = 9.2 x 10"^ h^^erg s'^), soft (0.5-2.0 
keV) X-ray sources at z ^ 0.2 on large scales, r — 5 — 100 
Mpc, f inding a correlati on length of tq — 7.4 ± 1.9 
h^^ Mpc. iGiin et all (pOOl measure the clustering of 



X-ray AGN at 2; = — 4 in both of the ~0.1 sq. degree 
Chandra Deep Fields (CDF), finding tq = 10.3 ± 1.7 h'^ 
Mpc and 7 — 1.33 ± 0.14 for 240 sources in the northern 
field and ro = 5.5 ± 0.6 h,-^ Mpc and 7 = 1.50 ± 0.12 
for 124 sources in the southern field. These error bars do 
not include cosmic variance, however, which is severe for 
the small CDF fields and most likely a ccounts for the dis - 
crepancy between the measurements. I Yang et"al] ()2006f ) 
measure the clustering of 233 spectroscopic sources at 
2 = 0.1 - 3 in the 0.4 sq. degree Chandra CLASXS 
well as 252 sources in the CDFN field. They 
find ro = 5.7 + 0.8/ - 1.5 h-^ Mpc in CLASXS and 
ro = 4.1 + 0.7/ - 1.1 h-^ Mpc in the CDFN field, 
where the errors are f rom bootst r ap res ampling of the 
data. More recently, 'Gilh et all ()2009[ ) find in the 2 
sq. degree XMM-COSMOS field that 349 X-ray AGN 
with 0.4 < z < 1.6 have a clustering scale length of 
ro = 5.2 ± 1.0 h-^ Mpc and 7 = 1.7 ± 0.2, where the 
error bars quoted include only Poisson errors. 

Ideally, one would like to measure the correlation am- 
plitude of AGN in a relatively narrow redshift range using 
spectroscopic redshifts, and compare with the clustering 
of galaxies in the same volume to understand the rela- 
tionship of AGN to their host galaxies and avoid much 
of the cosmic variance error. In this paper, we measure 
the clustering properties of 113 X-ray selected AGN in 
the redshift range 0.7 < z < 1.4 (a subset of the 463 
X-ray objects in this field with spectroscopic redshifts) 
using deep Chandra data in the Extended Groth Strip 
(EGS). For comparison, the only other X-ray field that 
has a large sample of spectroscopic redshifts and is com- 
parable in both size and depth to the EGS is the XMM- 
COSMOS field. The environments o f 53 X-ray selected 
AGN in the EGS were measured in iGeorgakakis et aTl 
(|2007f ) and compared to galaxies in the same volume; 
they found that X-ray AGN avoid under-dense regions 
and have environments similar to those of optically red 
galaxies. Here we measure the clustering of a larger sam- 
ple of AGN in the same field, as a function of scale, and 
compare to galaxies in the same volume. To measure 
the clustering amplitude of the AGN we do not use the 
auto-correlation function of the AGN, as it is very sen- 
sitive to spatially-varying selection functions (which is 
troublesome given that the Chandra point spread func- 
tion and telescope effective area vary across the field of 
view) and is subject to large shot noise due to the rel- 
atively small AGN sample size. Instead, we take ad- 
vantage of the large, uniformly selected DEEP2 galaxy 
survey in the EGS and measure the clustering of more 
abundant galaxies around the AGN using the two-point 
cross-correlation function and then divide by the corre- 
lation function of the galaxies themselves to determine 
the AGN clustering strength. We measure the cluster- 
ing of the AGN on both small (< 1 h^^ Mpc) and large 
(1 — 10 h~^ Mpc) scales and compare the results to the 
clustering of galaxies from the DEEP2 galaxy redshift 
survey as a function of color and luminosity at the same 
redshift, to understand which galaxies host these AGN 
and to constrain their host dark matter halo masses. By 
comparing these clustering measurements with those for 
optically-selected samples of luminous AGN, we can in- 
vestigate the relationship between accretion activity and 
large-scale structure, which strongly constrains theoret- 
ical models of the interplay between AGN and galaxy 
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Fig. 1. — Spatial distribution of DEEP2 galaxies (black squares) and X-ray AGN (colored diamonds) with 0.2 < z < 1.4 in the Extended 
Groth Strip field. X-ray AGN used for the clustering analysis performed here are shown in green; these sources are within the boundary 
of the DEEP2 survey in this field and cover the lower three quadrants of the strip, where the DEEP2 selection function is uniform (see 
Section 3.2 for details). 



evolution. 

The outline of the paper is as follows: §2 briefly de- 
scribes the Chandra and DEEP2 surveys in the EGS 
and presents additional spectroscopic follow-up of X-ray 
sources taken with the MMT. In §3 we discuss the optical 
properties of the X-ray AGN and define the various AGN 
and galaxy samples used here. §4 outlines the methods 
used to measure the cross-correlation function and infer 
the real-space auto-correlation function. Results for var- 
ious AGN and galaxy samples are shown in §5, as a func- 
tion of optical and X-ray properties. These results are 
discussed in §6, and we conclude in §7. To convert mea- 
sured redshifts to comoving distances along the line of 
sight, we assume a flat ACDM cosmology with ilm = 0.3 
and Qa = 0.7. We deflne h = -ffo/(100 km s"'^ Mpc-^) 
and quote correlation lengths, tq, in comoving Mpc. 

2. DATA 

2.1. Chandra Imaging 

X-ray AGN are identified using deep Chandra data 
(|Laird et al.ll2009| ) obtained as part o f the AEGIS survey 
of the Extended Groth Strip (EGS. [Pa^ et al.|[2007l ). 
This dataset consists of 8 Chandra ACIS-I pointings, 
each with a total integration time of about 200 ks split in 



at least 3 shorter exposures obtained at different epochs. 
The data reduction, source det ection and flux es tima- 
tion are described in detail by iLaird et "aTl (|2009l . In 
brief, after merging the individual observations into a 
single event file, images are constructed in four energy 
bands 0.5-7.0 keV (full), 0.5-2.0 keV (soft), 2.0-7.0 keV 
(hard) and 4.0-7.0 keV (ultra-hard). The count rates in 
the above energy intervals are converted to fluxes in the 
bands 0.5-10, 0.5-2, 2-10 and 5-10 keV, respectively. The 
limiting flux in each of these bands for a point source is 
estimated to be 3.5 x 10-l^ 1.1 x lO'^^, 8.2 x 10'^^, 
and 1.4 x 10~^^ ergs~^ cm~^, respectively. Sources are 
not required to be detected in both the soft and hard 
bands, though most are (for the main AGN sample used 
here, described in Section 3 below, 29% of the sources 
are detected in the soft band only and 9% are detected 
in the hard band only). The X-ray catalog comprises a 
total of 1325 unique sources over 0.63 deg^ to a Poisson 
detection probability threshold of 4 x 10~^. We identify 
optical counterparts of the Chandra X-ray sources using 
the maximum like lihood (LR) method as described in 
ILaird et all (|2009f ). A total of 895 sources have optical 
counterparts to Rab = 24.1 mag {LB. > 0.5). 
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Fig. 2. — Left: Optical Rab magnitude distribution for X-ray sources targeted (dotted line) for MMT/Hectospec spectroscopic observa- 
tions. The solid line shows objects which yielded a high-confidence spectroscopic redshift. Right: Optical apparent B — R colors for X-ray 
sources targeted (dotted line) and for those with a high-confidence redshift (solid line). 
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Fig. 3. — Left: Redshift distribution of all 463 Chandra X-ray sources in AEGIS with spectroscopic redshifts (solid line) and DEEP2 
galaxies in the EGS (dashed line). The dotted vertical lines show the redshift range used here to measure the clustering properties of X-ray 
AGN and compare with galaxies. Right: The redshift distribution of X-ray sources where the redshift was measured in Keck/DEIMOS 
spectra (solid line, 233 sources) or MMT/Hectospec spectra (dashed line, 230 sources). 

2.2. DEEP2 Spectroscopy 

Spectroscopic redshifts were obtained for 463 X-ray 
sources in AEGIS, 233 (50%) of w hich are from the 
DEEP2 Galaxy Redshift Survey (Da vis et al.ll200ai2003f) 
and 230 (50%) of which were obtained with follow-up 
spectroscopic observations at the MMT, as described in 
§2.3. The DEEP2 survey is a com pleted project using the 
DEIMOS spectrograph (|Faber e t al. 2003) on the 10m 
Keck II telescope to survey optically-selected galaxies at 
2; ~ 1 in a comoving volume of approximately 5x10® h^^ 
Mpc^. Using ~ 1 hour exposure times, the survey has 
measured high-confidence redshifts for > 30, 000 galax- 
ies, most in the redshift range 0.7 < z < 1.5, to a limiting 
magnitude of Rab = 24.1. Roughly 65% of the galax- 
ies to i?AB = 24.1 were targeted for spectroscopy, and 
of those targeted, successful redshifts were measured for 
-65%. 

One of the DEEP2 fields is the EGS, where the sur- 
vey has measured > 11,000 high-confidence redshifts 



(jPavis et al.ll2007f l: unlike other DEEP2 fields, objects 
with colors indicating that they have z < 0.7 are still ob- 
served in the EGS. The DEEP2 spectra have high resolu- 
tion (i? ~ 5, 000), and rms redshift errors (as determined 
from repeated observations) are < 35 km s^^. After re- 
stricting to 0.7 < z < 1.4 and the area covered by the 
AGN samples described below, the EGS DEEP2 sample 
contains a total of 4669 sources. 

2.3. MMT Spectroscopy 

To obtain redshifts for more X-ray AGN than were 
targeted in the DEEP2 survey (which began before 
the X-ray observations were taken, so X-ray sources 
were not given higher observing priority), we used the 
MMT/Hectospec fiber spectrograph for follow-up spec- 
troscopy of optical counterparts to X-ray sources. We 
observed 5 Hectospec configurations in queue mode in 
May 2007, July 2007, and May 2008. Total integration 
times were 1.5-2 hours per configuration. The wave- 
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length coverage was ~ 4500 — 9000 A at 6 A resolu- 
tion. Observing conditions were fair, though not ideal; 
one night was affected by moonlight, and the seeing was 
often > 1". Additional targets were observed as filler ob- 
jects on configurations taken by other AEGIS team mem- 
bers (C. Willmer and P. Barmby, priv. comm.) during 
the same observing seasons, with similar spectrograph 
setups. In total, we targeted optical counterparts for 449 
X-ray sources. 

The data were reduced using the HSRED IDL re- 
duction pipeline^^. Redshifts were measured using a 
^^-minimization between the observed spectrum and 
emission-line and absorption-line galaxy and AGN tem- 
plates, and were confirmed by eye. We obtained high- 
confidence redshifts for 230 sources, or 51% of the tar- 
geted sample. 

The spatial distribution of the X-ray AGN with red- 
shifts between 0.2 < z < 1.4 is shown in Figure [1] (col- 
ored diamonds), along with the distribution of DEEP2 
galaxies in the same redshift range (black squares). As 
discussed below in Section 3.2, galaxies and AGN in the 
upper quadrant are not used for the clustering measure- 
ments performed here; the green squares in Figure [1] show 
the X-ray AGN used to measure the cross-correlation 
function. 

The i?AB magnitudes and B ~ R colors of the targeted 
sources and the subsample with successful redshift mea- 
surements are shown in Figure [H Our success rate de- 
clines at i?AB > 22.5, as expected, as the signal-to-noise 
in the 2-hour spectra was not high enough to measure 
a robust redshift. A K-S test of the B — R color distri- 
bution for the targeted sources and the subsample with 
successful redshifts gives a 93% probability that they are 
drawn from the same parent population. The success 
rate is lowest at B — R ^ 1, but is comparable for red 
and blue galaxies. 

Figure [3] shows the redshift distribution for all the X- 
ray sources with spectroscopic redshifts in the EGS (left), 
and split according to the instrument used to obtain 
spectroscopy (right). The left panel also shows the red- 
shift distribution of DEEP2 galaxies in the EGS (dotted 
line), scaled down by a factor of 35. The redshift distri- 
bution for X-ray sources observed with Keck/DEIMOS 
(right panel) shows a peak dX z ^ 0.7 and a small tail 
of objects at z > 1.5, reflecting both the underlying red- 
shift distribution of the i?AB < 24.1 DEEP2 sources and 
the high resolution of the spectra, which results in a nar- 
rower spectral range of ~ 6500 - 9000 A, limiting the 
redshifts where multiple features will be observed. The 
MMT/Hectospec redshift distribution is flatter and in- 
cludes more objects at z > 1.5, due to the broader wave- 
length coverage of the spectra. 

K-corrections, absolute Mb magnitudes and restframe 
(U — B) colors ha ve been derived as described in 
I Willmer et aI](|2006D . Absolute magnitudes given in this 
paper are in the AB system and are Mb — 5 \og{h) with 
/i = 1, which for the remainder of the paper we denote 
as Mb- X-ray luminosities are derived using h = 0.7. 

3. AGN PROPERTIES AND SAMPLES 
3.1. Optical Properties of X-ray AGN 

See |http: / /www.astro.princeton.edu/~rcool/hsred| 
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Fig. 4. — Absolute magnitude Mb as a function of redshift for 
DEEP2 galaxies in AEGIS (black crosses) and X-ray AGN (red 
diamonds). The AGN typically reside in optically bright galax- 
ies, though they are occasionally found in fainter galaxies. AGN 
subsamples as a function of Mg are drawn to test the luminosity- 
dependence of AGN clustering. 

Figure m shows the absolute magnitude, Mb, for the X- 
ray AGN host galaxies (red diamonds) and all DEEP2 
galaxies (black crosses) as a function of redshift. The 
AGN tend to reside in optically bright galaxies, along 
the upper portion of the distribution, at all redshifts be- 
tween 2: = 0.2 — 1.4. Figure [5] shows the color-magnitude 
diagram for the X-ray AGN (red diamonds) and DEEP2 
galaxies (black crosses) in the redshift range z — 0.2 — 1.4 
(le ft) and z = .7 — 1.4 (right). As originally shown 
in iNandra et al.l (|2007f) and update d here with a larger 
sample (see also lBarger et al.l (|2003f )l. X-ray AGN in the 
200-ks depth Chandra AEGIS data predominantly have 
host galaxies that are either on the red sequence or the 
bright /red part of the blue cloud, implying that they re- 
side in massive host galaxies (see also Silverman et al. 
2008). X-ray AGN that are very bright and blue, with 
Mb < —21.3 and {U ~B) < 0.3, in the lower right corner 
of the figure, we define to be quasars, where the optical 
light is dominated by the AGN and not the host galaxy. 
ACS imaging of these sources in the EGS shows a bright 
blue point source with no host galaxy visible (Pierce et 
al. 2007, Georgakakis et al. in preparation) . . 14% of the 
X-ray AGN in our sample at 0.7 < z < 1.4 are quasars, 
according to this definition. We note that contamination 
of the X-ray AGN sample by normal star-forming galax- 
ies witho ut AGN in this sample is expected to be less 
than 1% (|Montero-Dorta et al.|[2008l ). 

In Figure [6] we compare the optical {U — B) color 
with the hardness ratio of the X-ray source, measured 
as HR={H- S)/{H + 5), where = 2 - 7 keV counts 
and S = 0.5 — 2 keV counts. The hardness ratio esti- 
mates used here are determined through Bayesian meth- 
ods and do not require detection in both the hard and soft 
bands; the actual source counts and background counts 
in both bands are used even if the object is not signifi- 
cantly detected. Therefore none of the X-ray AGN have 
HR = — 1 or HR = +1. The dashed line roughly marks 
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Fig. 5. — Restframe optical color-magnitude ({U — B) versus Mb) diagram for DEEP2 galaxies (black crosses) and X-ray AGN (red 
diamonds) in AEGIS for the redshift ranges z = 0.2 — 1.4 (left) and z = 0.7 — 1.4 (right). The dashed line shows the definition of the 
minimum in the restframe {U — B) bimodality, used to define red and blue galaxies. X-ray AGN are predominantly found in red sequence 
galaxies, 'green' galaxies in the dip of the color bimodality, and the brightest/reddest of the blue cloud galaxies. The optical light for a 
small percentage (^ 14%) of the X-ray sources is dominated by the AGN itself and not the host galaxy; these sources are seen in the lower 
right corner and are optically very bright and blue. 
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Fig. 6. — Left: Optical (U — B) color as a function of hardness ratio (see §3.1 for details) for all X-ray AGN in the redshift range 
z = 0.2 — 1.4. The dashed line roughly indicates the separation of red and blue galaxies, above and below the line respectively. The 
optically red and 'green' transition galaxies show a range of hardness ratios, while the blue galaxies are on average softer sources. Right: 
Optical absolute magnitude Mb as a function of hardness ratio for all X-ray AGN in the redshift range 2: = 0.2 — 1.4. A small but significant 
correlation is found between optical magnitude and hardness ratio. 



the minimum of the color bimodahty, with blue galaxies 
at {U - B) < 1.0 and red galaxies at {U - B) > 1.0. 
We find a correlation between hardness ratio and op- 
tical obscuration in t hat the hardest X-r ay sources are 
all optically red fe.g.. lMainieri et al]|2002[ ). This may ei- 
ther be the result of obscuration in the host galaxy itself, 
which reddens the galaxy as well as obscures the AGN, 
or it may be that obscured AGN only reside in galax- 
ies with older stellar population s (|Nandra et al.l l2007t 
iRovilos &: Georgantopoulosll2007D . We find that the red 
and 'green' galaxies have a wide range of hardness ratios, 
while the blue galaxies are softer sources. A small but 
significant correlation is found between absolute mag- 
nitude Mb and hardness ratio, which have a Spear- 
man's rank correlation coefficient oi p = 0.14 at 97% 
significance. This correlation persists after excluding the 



brightest sources with Mb < —22, indicating that it is 
not likely to be dominated by an increased contribution 
from unobscured AGN. 

3.2. AGN samples 

From the full AGN and galaxy datasets we define sev- 
eral samples used to measure the AGN-galaxy cross- 
correlation function. We limit most samples to the red- 
shift range 0.7 < z < 1.4 so as to match the redshift range 
of the other DEEP2 fields where we have performed ex- 
tensive clustering an alyses of galaxies as a function of 
color and luminosity (|Coil et al.ll2006al . |2008[ ). as well as 
the measur ement of the q uasar-galaxy cross-correlation 
function (Coil et al1l2007f ). with which we compare our 
results here. We do not include data from the upper 
quadrant of the EGS field (pointing 14 in the DEEP2 



7 



45 



44 



43 



42 



41 



40 








* 










V"'*'"*"^" 

\ : ' ■ 


















1,1.1, 





S! 0.0 



o 

I 



0.4 



0.6 0.8 

Redshift 



1.0 



1.2 




0.6 0.8 

Redshift 



Fig. 7. — Left: X-ray luminosity as a function of redshift for all X-ray AGN in the redshift range z = 0.2 — 1.4. The 'AGN high Lx' and 
'AGN low Lx' samples are defined by the dashed and solid lines, respectively. Right: Hardness ratio of the X-ray AGN as a function of 
redshift. The 'AGN hard' and 'AGN soft' samples are defined as being above or below the solid line, corresponding to Njj = 10^^ cm~^ 
for r = 1.9. 



TABLE 1 
AGN AND Galaxy Samples 



Sample No. of 
objects 


z median 
range z 


median 
Mb 


selection^ 


AGN main 


113 


0.7-1.4 


0.90 


-20.98 


does not include quasars 


AGN host bright 


54 


0.7-1.4 


1.00 


-21.33 


Mb < -21 


AGN host faint 


59 


0.7-1.4 


0.84 


-20.45 


Mb > -21 


AGN host red 


50 


0.7-1.4 


0.90 


-20.96 


red iU — B) host galaxy 


AGN host blue 


63 


0.7-1.4 


0.91 


-20.98 


blue (U — B) host galaxy 


AGN low Lx 


91 


0.4-1.2 


0.74 


-20.40 


42 < logLx < 43 erg s-l 


AGN high Lx 


83 


0.4-1.2 


0.84 


-21.01 


log Lx > 43 erg 


AGN hard 


140 


0.2-1.4 


0.75 


-20.62 


Nh > lO^^cm.-^ 


AGN soft 


111 


0.2-1.4 


0.67 


-20.60 


Nh < 1022cm-2 


All galaxies 


4669 


0.7-1.4 


0.92 


-20.16 




Red galaxies 


757 


0.7-1.4 


0.85 


-20.66 


red {U - B) 


Blue galaxies 


3913 


0.7-1.4 


0.94 


-20.06 


blue {U - B) 


Matched galaxies 


504 


0.7-1.4 


0.98 


-20.97 


Red matched galaxies 


226 


0.7-1.4 


0.92 


-20.92 




Blue matched galaxies 


278 


0.7-1.4 


1.03 


-20.99 





" See Section 3 for details on each AGN sample and Section 5 for details on each 
sample. 



ilaxy 



naming convention), where the DEEP2 target selection 
function is shghtly different from the rest of the survey, 
due to shallower BRI photometric imaging l|Coil et al.l 
[2004a) used in the target selection. The total area cov- 
ered by the data used here is ~ 0.4 deg^. 

Details of each AGN sample are given in Table [T] The 
'AGN main' sample includes all non-quasar X-ray AGN 
over the redshift range z — 0.7 — 1.4, a total of 113 
sources; as seen in the left panel of Figure [7] these AGN 
have log Lx > 42 erg s~^, where Lx is the restframe 
2-10 keV luminosity calculated from the observed 2-10 
keV flux, assuming a photon index F = 1.9. 

The main AGN sample has a mean log Lx = 42.4 erg 
and a median log Lx = 42.8 erg The 'AGN 

host bright' and 'AGN host faint' samples are defined 
by their optical Ms magnitude; the samples are split at 
the median Mb magnitude of the AGN main sample. 
The 'AGN host red' and 'AGN host blue' samples are 
defi ned by the optical co lor of the host galaxy. Follow- 
ing iWillmer et all (l2nn6h . we define red and blue galax- 



ies using the observed color bimodality in DEEP2 with 
the following tilted cut in color-magnitude space (in AB 
magnitudes): 



{U -B) = -0.032(Ms + 21.62) + 1.035, 



(1) 



as shown in Figure [5l We do not allow this color cut to 
evolve with redshift within the redshift range used here. 

For the AGN samples defined by optical color or mag- 
nitude we exclude the quasars in the lower right portion 
of Figure [3 as the optical color and magnitude of these 
objects is not reflective of the host galaxy but is from 
the AGN itself. We define quasars here as sources with 
Mb < —21.3 and {U — B) < 0.3; there are 19 quasars in 
our sample with 0.7 < z < 1.4 (12 of which are shown in 
the right panel of Figure \5\ the rest are outside the plot 
range), out of 132 X-ray sources. None of our results 
change if these sources are included. For the 'AGN host 
bright/faint' and 'AGN host red/blue' samples, the opti- 
cal magnitude and color are properties of the host galaxy, 
not the AGN itself. The sample of spectroscopically- 
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identified broad-line AGN is too small to use for robustly 
measure clustering properties; a larger sample is needed. 

We also divide the full AGN sample by X-ray lumi- 
nosity, using the median value of log = 43 erg s~^, 
such that the 'AGN low L^' sample is defined as having 
42 <log Lx < 43 erg (which results in a roughly 
volume-limited sample), while the 'AGN high Lx' sam- 
ple has log Lx > 43 erg s~^, as shown in the left panel 
of Figure [71 Since here we are concerned with compar- 
ing the Lx samples to each other and not to galaxies 
in the same redshift range, we are free to use a differ- 
ent redshift range from that used in the DEEP2 galaxy 
clustering studies. For the Lx-selected samples we use a 
lower redshift limit of z = 0.4, which increases the sample 
size but is not so low as to cause the two samples to have 
significantly different mean redshifts (see Table 1). To 
keep the 'AGN low Lx' sample roughly volume-limited 
we use an upper redshift limit of z = 1.2 (see Figure 
[7]), which ensures that any difference seen in their clus- 
tering properties should be dominated by differences in 
luminosity and not redshift. The 'AGN hard' and 'AGN 
soft' samples are defined by the X-ray color, or hardness 
ratio, as shown in the right panel of Figure [71 The solid 
line corresponds to Nh = lO^^cm"^ for F = 1.9 and de- 
fines the 'AGN hard' and 'AGN soft' samples. For these 
samples we use the full redshift range of z = 0.2 — 1.4. 

Several galaxy samples are used in this paper; details 
for each are given in Table 1. All cross-correlation func- 
tions presented here use the full DEEP2 galaxy sample 
with 0.7 < z < 1.4 ('AH galaxies' in Table 1) to cross- 
correlate with either an AGN or a galaxy sample. This 
full galaxy sample is separated into red and blue galaxy 
samples based on restframe {U — B) color, with the color 
cut given above. Matched galaxy samples are discussed 
in §5.3; these samples have the same Mb and ([/ — B) 
distributions as various comparison AGN samples. 

4. MEASURING THE CROSS-CORRELATION FUNCTION 

To measure the clustering of the X-ray AGN sample, 
we use the cross-correlation of AGN with galaxies in the 
same volume, and the auto-correlation function of the 
galaxies themselves, to infer the auto-correlation function 
of the AGN sample. 

The two-point auto-correlation function ^(r) is defined 
as a measure of the excess probability above that for an 
unclustered distribution of finding an object in a volume 
element dV at a separation r from another randomly 
chosen object, 

dP ^ n[l + i{r)]dV, (2) 

where n is the mean n umber density of the object in 
question ([Peebles! 1 19801) . The cross-correlation function 
is the excess probability above Poisson of finding an ob- 
ject from a given sample at a separation r from a random 
object in another sample. Here we measure the cross- 
correlation between AGN and galaxies: 

dPiG\A)^nG[l + UGir)]dV, (3) 

which is the probability of finding a galaxy (G) in a vol- 
ume element dV at a separation r from an AGN (A), 
where is the number density of galaxies. 

To estimate the cross-correlation function between our 
AGN and galaxy samples, we measure the observed num- 
ber of galaxies around each AGN as a function of dis- 
tance, divided by the expected number of galaxies for a 



random distribution. We use the estimator 



where AG are AGN-galaxy pairs and AR are AGN- 
random pairs at a given separation, where the pair counts 
have been normalized by and nn, respectively, the 
mean number densities in the full galaxy and random 
catalogs. This estimator is preferred here as it does not 
require knowledge of the AGN selection function, only 
the galaxy selection function, which is well-quantified. 
To calculate ^(r) we create a random catalog with the 
same redshift distribution as the DEEP2 galaxies used 
for the cross-correlation measurement and the same sky 
coverage as the DEEP2 galaxies in the EGS, applying the 
two-dimensional variation in the DEEP2 target selection 
and redshift success rate over the plane of the sky. We 
also mask the regions of the random catalog where the 
photometric data are affected by saturated stars or CCD 
defects. Details on the random ca talog construction are 
given in §3.1 of iCoil et al.l (|2004bD . 

We wish to measure the real-space correlation function, 
^(r). However, peculiar velocities distort the positions 
of objects in redshift space along the line of sight. To 
capture this effect, we measure ^ as a function of two 
coordinates, perpendicular to (rp) and along (tt) the line 
of sight. In applying the above estimator to galaxies, 
pair counts are computed over a two-dimensional grid of 
separations to estimate ^{vp, tt). To recover ^(r), ^(rp, tt) 
is integrated along the tt direction and projected onto 
the Tp axis. As redshift-space distortions affect only the 
line-of-sight component of ^(rp, tt), integrating over the tt 
direction leads to a statistic Wp{rp), which is i ndependent 
of redshift-space distortions. Following iDavis fc Peeblea 
11983), 

Wp{rp) = 2 / dTT ^(rp, tt) - 2 / dy ^[(r^ + 
Jq Jo 

(5) 

where y is the real-space separation along the line of 
sight. Here, as in other DEEP2 clustering studies, we 
integrate to a maximum separation in the vr direction of 
20 Mpc, as the signal to noise degrades quickly for 
larger separations where ^ becomes small. Systematic 
e ffects due t o the use of slitmasks are small (see §3.4 
of iCoil et al.l (|2008[ ) for details). We do not correct for 
these effects here, as we are primarily concerned with 
the relative clustering of AGN samples compared to each 
other and to the clustering of galaxy samples, and to first 
order these effects will cancel. 

Error bars on Wplrp) are estimated from jackknife sam- 
pling of the data, where the EGS is divided into eight 
regions of equal area and Wp(rp) is calculated for eight 
jackknife samples, removing one area at a time. To check 
the validity of these jackknife error bars, we compare 
the Wp{rp) errors for the 'main AGN sample' with er- 
rors that include both Poisson and cosmic variance er- 
rors estimated using DEEP2 mock catalo gs derived from 
the M illenium Run simulation fsee .Kitzbichler k, White! 
[20071 for further details). Using a total of 24 indepen- 
dent mock catalogs with the same geometry and selection 
function as the DEEP2 data in the EGS, we select 113 
random red galaxies with 0.7 < z < 1.4 to be X-ray 
AGN proxies, since red DEEP 2 galaxies have a similar 
clustering amplitude as the X-ray AGN, as shown below 
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TABLE 2 
Relative Bias Results 



Samples Relative bias Relative bias 

0.1 < rp < 8 Mpc 1 < Tp < 8 Mpc 



AGN host bright /faint 

AGN host red/blue 

AGN high Lx/low Lx 

AGN hard/soft 

AGN/Red galaxies 

AGN/Blue galaxies 

AGN/Matched galaxies 

AGN host red/Red galaxies 

AGN host red/Red matched galaxies 

AGN host blue/Blue galaxies 

AGN host blue/Blue matched galaxies 



1.52 ± 0.32 


1.20 ±0.31 


1.91 ± 0.38 


1.79 ±0.27 


1.05 ± 0.18 


0.95 ± 0.20 


1.22 ± 0.12 


1.27 ±0.17 


0.97 ±0.07 


0.94 ±0.08 


1.61 ± 0.11 


1.48 ±0.12 


1.28 ± 0.10 


1.15 ±0.09 


1.20 ±0.09 


1.13 ±0.09 


1.33 ± 0.15 


1.19 ±0.12 


1.27 ±0.15 


1.22 ±0.13 


1.24 ±0.16 


1.09 ±0.16 



in Section 5.3. We calculate the cross-correlation of these 
AGN proxy galaxies with all mock DEEP2-like galaxies 
in each mock catalog, and then calculate the variance 
across the 24 catalogs. The error on Wp{rp) in the mock 
catalogs is less than the error estimated from jackknife 
sampling of the data itself: on scales rp=l-8 Mpc in 
the mock catalogs the fractional error on Wp{rp) is 8%, 
while in the data it is 10%. To be conservative we there- 
fore use the error estimated from the jackknife sam ples 
of the data. As previously noted (jCoil et al.n 2007\ the 
Millenium Run mock catalogs appear to exhibit smaller 
variance than the real Universe. 

If ^(r) is modeled as a power law, ^(r) = (r/rg)"''', 
then rp and 7 can be readily extracted from the projected 
correlation function, Wp{rp), using an analytic solution to 
Equation [5] 



Wp{rp) 



ro 



r(^)r(V) 
r(i) 



(6) 



where F is the usual gamma function. A power-law 
fit to Wpij-p) will then recover vq and 7 for the real- 
space correlation function, ^(r). In practice, however, 
we can not measure S.irpjir) accurately to infinite sep- 
arations as assumed in Equation [5l and here we inte- 
grate Wp{rp) to TTijjax = 20 h^^ Mpc. To recover the 
correlation length and slope wc then estimate Wp{rp) in- 
tegrated to n,nax — 20 h^^ Mpc for a grid of rg and 
7 values and minimize the difference betwee n the data 
and model. Full details are given in §4.2 of iCoil et al.l 
([2008). We do not use a full covariance matrix in the 
fits. Errors on tq and 7 are propagated from the jack- 
knife errors on Wp{rp); if instead we measure rp and 7 
in the separate jackknife samples and use the variance 
among the various fits, which does take into account the 
covariance between the Vp bins, the resulting errors are 
somewhat smaller. To be conservative we use the larger 
errors derived from the Wp{rp) errors. 

5. RESULTS 

5.1. Clustering oj Main AGN Sample 

The projected cross-correlation of the 'AGN main' 
sample with all DEEP2 galaxies in the same redshift 
range, Wp{rp), is shown with a dashed red line in Fig- 
ure [51 As discussed above, the Wp{rp) error bars are de- 
rived from jackknife samples of the data and are slightly 
larger than errors inferred from cosmological simulations. 
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Fig. 8. — The projected AGN-galaxy cross-correlation func- 
tion Wp{rp) as a function of scale for the non-quasar X-ray AGN 
main sample (red dashed line) and the full DEEP2 galaxy popula- 
tion in the EGS with 0.7 < z < 1.4. The galaxy auto-correlation 
function in the same volume is shown (blue dotted line), as well 
as the inferred X-ray AGN auto-correlation function (black solid 
line). Errors are computed from jackknife samples of the data and 
are comparable to errors found from mock catalogs (see text for 
details). 

The projected galaxy auto-correlation function is also 
shown (blue dashed line) for the full galaxy sample with 
0.7 < z < 1.4. Assuming a linear bias, 



Wp{AG) = [wpiAA) X Wp{GG)]^^^, 



(7) 



where Wp{AG) is the AGN-galaxy cross-correlation func- 
tion, and Wp{AA) and wp{GG) are the AGN and 
galaxy auto-correlation functions. From measurements 
of 'Wp{AG) and wp{GG), as shown in Figure [51 the X- 
ray AGN auto-correlation function wp{AA) can then be 
inferred; this is shown in Figure [8] as a solid black line, 
where the errors are derived from jackknife resampling 
of wp{AG)^/wp{GG). 
Fitting Wp{AA) as a power law, we find tq = 5.95±0.90 
Mpc and 7 = 1.66 ± 0.22. To compute the bias 
of the X-ray AGN auto-correlation function compared 
to that expected for dark matter, we estimate Wp(rp) 
for dark matter particles at the mean redshift of the 
main AGN s ample , using the publicly-available code of 
[Smith et al.l (|2003[ ). where we integrate the dark mat- 
ter C(rp,7r) to TTmax = 20 Mpc. We estimate the 
linear galaxy bias from the ratio of the quantities — 
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Fig. 9. — The projected AGN-galaxy cross-correlation functions for AGN samples split according to optical Mg magnitude (upper left), 
optical {U — B) color (upper right), X-ray luminosity (lower left), and X-ray color or hardness ratio (lower right). Details of each sample 
are given in Table 1. No significant difference in clustering is detected for the samples selected on optical luminosity, X-ray luminosity or 
X-ray color. There is a correlation with optical color: the X-ray AGN with red host galaxies are significantly more clustered than those in 
blue host galaxies. 



[wj,]gai/[wp]darkmatter, and average over scales 1 < < 8 
h^^ Mpc, to find a mean bias oi b = 1.85 ± 0.28. 

From the bias of the main AGN sample we can 
infer the minimum dark matter mass for halos that 
host these AGN. A s sho wn in the Appendix of 
IZheng. Coil, fc Zehavil (|2007l ). the dark matter mass in- 
ferred assuming that all objects, whether galaxies or 
AGN, are 'central' objects in their dark matter halos 
(ie., assuming one galaxy or AGN per halo) can lead to 
an overestimate of the minimum halo mass. Here we 
assume that the AGN population has the same satel- 
lite fraction, or number of additional galaxies per halo, 
as found in the DEEP2 galaxy population using HOD 
modeling of the co rrelation function, where the satellite 
fraction is ~15% (|Zheng. Coil, fc Zehavil [20071 ) . Using 
this correction for the satellite fraction changes the rel- 
ative weights given to halos of a different mass when 
computing the large-scale bias for different halo masses, 
and the minimum mass inferred is slightly lower than as- 
suming one AGN per halo. Here we find that a± z — 0.94 
the minimum dark matter halo mass of the non-quasar 



X-ray AGN is M™„, = 5 (-1-5/ - 3) x W'^h-^M^,, sim ilar 
to red galaxies at the same redshift (|Coil et al.ll20?)8l ). 

5.2. Clustering of Other AGN Samples 

To investigate the dependencies of AGN clustering 
properties, we divide the main AGN sample into sam- 
ples based on optical luminosity, optical color. X-ray 
luminosity and hardness ratio. Details of the selection 
for each sample are given in Section [3] and in Table [TJ 
Figure [9] shows the projected cross-correlation with the 
DEEP2 galaxies, Wp{rp), for each of these samples. To 
quantify differences in clustering amplitude, we calcu- 
late the relative bias between these samples. The rela- 
tive bias is defined here as the ratio of Wp{rp) for two 
samples. Note that this is the usual meaning of relative 
bias {b^^i = ^1/^2), as here we are using cross-correlation 
functions with the identical galaxy tracer samples and 
not auto-correlation functions. The error bars on the 
relative bias measurement are not estimated from the 
errors on Wp{rp) directly, as that would not include co- 
variance between the bins and would include cosmic vari- 
ance errors, which cancel to first order when comparing 
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Fig. 10. — The projected AGN-galaxy cross-correlation function 
Wp{rp) as a function of scale for the AGN main sample (black 
solid line) and the full DEEP2 galaxy sample with 0.7 < z < 1.4. 
The cross-correlations of red and blue DEEP2 galaxies with the 
full galaxy sample in the same volume are also shown with dashed 
red and dotted blue lines, respectively. The X-ray AGN are seen 
to have very similar clustering properties to red galaxies and are 
more clustered than blue galaxies at 2 ~ 1. 

the clustering of two populations in the same volume. 
The relative bias errors are therefore estimated from the 
variance of the relative bias measured across the jack- 
knife samples. 

As the relative bias can be scale-dependent, we calcu- 
late the unweighted mean relative bias over two scales, 
rp = 0.1-8 h^^ Mpc ('all scales') and = 1 - 8 /i"^ 
Mpc ('large scales'), and also measure the variance of 
the relative bias on the same scales. Relative bias re- 
sults for various samples are given in Table [21 We find 
no significant difference in the clustering of AGN when 
selected by optical luminosity, X-ray luminosity, or hard- 
ness ratio, within the ranges sampled here. We find that 
the AGN in red host galaxies are significantly more clus- 
tered than those in blue host galaxies, at the 2.9a level 
on large scales, with a relative bias of brei = 1.79 ± 0.27, 
and at the 2Aa level on all scales, with a relative bias of 
brei = 1.91±0.38. This relative bias is consistent with the 
relative bias found between red a nd blue DE EP2 galaxies 
of brei = 1.44 ± 0.07 measured in lCoil et all ([2008). 

We do not find a significant difference in the clustering 
of hard X-ray sources compared to soft X-ray sources. 
This lack of a detectable correlation between hardness 
ratio and clustering amplitude may be surprising, given 
that we observe correlations both between host {U — B) 
color and hardness ratio (in that the hardest sources are 
in optically red host galaxies, see Figure [6]) and between 
host [U — B) color and clustering amplitude, both us- 
ing smaller samples than for the hardness ratio-clustering 
comparison. One reason for the lack of a significant cor- 
relation between hardness ratio and clustering amplitude 
may be the large range of hardness ratios among the red 
AGN host galaxies (see Figure [H]), so that while the hard 
sample is dominated by red galaxies, the soft sample has 
both red and blue galaxies, mitigating the effect of the 
hardness-color and color- clustering correlations. Our re- 
sults clearly show that AGN host color is more correlated 



Fig. 11. — The projected AGN-galaxy cross-correlation func- 
tion Wp(rp) for the AGN main sample compared to the cross- 
correlation function of a sample of galaxies with the same optical 
color and magnitude distribution, regardless of whether they host 
AGN. Those galaxies with X-ray AGN are found to cluster more 
than the full galaxy population with the same color and luminosity 
distribution, which implies that they reside in more massive dark 
matter halos. 

with dark matter halo mass than hardness ratio is: i.e., 
star formation histories are more closely related to local 
environment than AGN X-ray hardness is. 

5.3. Clustering of AGN Compared to Galaxies 

In Figure [TO] we compare the main AGN sample 
with the clustering of red and blue DEEP 2 galaxies 
in the same volume. Here we have computed the 
cross-correlation of red or blue DEEP2 galaxies with 
all DEEP2 galaxies, similar to the AGN-galaxy cross- 
correlation measurement. We use all red and blue 
DEEP2 galaxies, without an absolute magnitude limit, 
to ensure that the redshift distribution of the red and 
blue galaxies is similar to that of the X-ray AGN. The 
median absolute magnitudes of the galaxy samples are 
Mb = -20.1 for the blue galaxies and Mb = -20.7 for 
the red galaxies. If we do impose a limit of Mb < —20 
for the blue galaxies it does not change the result; the rel- 
ative bias is within la of the value using all blue DEEP2 
galaxies. 

The X-ray AGN are seen to have a similar clustering 
amplitude as red DEEP2 galaxies and are more clustered 
than blue DEEP2 galaxies. The relative bias of the main 
AGN sample to red galaxies is brei = 0.94 ± 0.08 on 
large scales and brei = 0.97 ± 0.07 on all scales, while 
the relative bias to blue galaxies is brei = 1-48 ± 0.12 on 
large scales and brei = 1.61±0.11 on all scales. The X-ray 
AGN are significantly more clustered than blue galaxies 
at z ~ 1. 

The slope of the X-ray AGN correlation function is 
more similar to t hat of blue DEE P2 galaxies than red 
DEEP2 galaxies (jCoil et al.ll2008l) . but within 2a it is 
consistent with either. As seen in Figure [3 the X-ray 
AGN correlation function may have a steeper slope on 
small scales, rp < 0.5 Mpc, compared to all galaxies, 
but the error bars on the X-ray AGN Wp(rp) are too large 
to detect a significant difference. 

While the main AGN sample has the same clustering 
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Fig. 12. — The clustering of AGN host galaxies compared to all galaxies of a given color. Left: The projected AGN-galaxy cross- 
correlation function Wp{rp) as a function of scale for AGN in red host galaxies (solid line) compared to all red galaxies (dashed line) and red 
galaxies with the same color and magnitude distribution (dotted line) at 0.7 < z < 1.4. Right: The projected AGN-galaxy cross-correlation 
for AGN in blue host galaxies (solid line) compared to all blue galaxies (dashed line) and blue galaxies with the same color and magnitude 
distribution (dotted line) at 0.7 < z < 1.4. Galaxies which host AGN are found to be somewhat more clustered than galaxies of the same 
color and magnitude, especially on scales rp < 1 h^^ Mpc. 



amplitude as red galaxies at z ~ 1, only half (42%) of 
X-ray AGN are in red host galaxies. This implies that 
AGN in red galaxies may be more clustered than red 
galaxies as a whole and/or that AGN in blue galaxies 
may be more clustered than blue galaxies as a whole. It 
also implies that galaxies of a given color and luminosity 
distribution that host X-ray AGN may be more clustered 
than those without AGN. 

We can test this by comparing the clustering of the 
main AGN sample with that of DEEP2 galaxies with the 
same color and magnitude distribution as the AGN host 
sample; the results are shown in Figure 1111 To create 
the 'matched galaxy sample' to compare with the main 
AGN sample, we first measure the joint color-magnitude 
distribution of the main AGN sample in {U — B) and 
Mb space, in bins of 0.1 in {U — B) color and 0.3 in 
Mb magnitude. We then normalize this distribution to 
have a maximum value of 1, making it effectively a rela- 
tive weight map in color-magnitude space. We then se- 
lect all DEEP2 galaxies in the same volume as the main 
AGN sample, and for each bin in color-magnitude space 
we randomly sample the galaxy population (without re- 
placement) such that the fraction of the final sample in 
that bin is equal to the weight in the main AGN sample. 
We have tested that creating a matched sample with re- 
placement from the parent population does not change 
the results within the \a error bars. The matched sam- 
ple therefore has a nearly identical joint distribution in 
color-magnitude space as the main AGN sample. The re- 
sulting cross-correlation of this 'matched galaxy sample' 
with all DEEP2 galaxies in the same volume is shown as 
the dashed line in FigurefTTI The AGN are seen to cluster 
more than the matched galaxy sample at the 2.8(t level 
when measured on scales 0.1 < rp < 8 Mpc, with a 
relative bias of Kei = 1-28 ± 0.10, and Kei = 1-15 ± 0.09 
on scales 1 < < 8 h^^ Mpc. Galaxies of a given opti- 
cal color and luminosity which host an AGN are therefore 
more clustered than their quiescent counterparts. 

The higher clustering amplitude found for AGN in red 



galaxies compared to AGN in blue galaxies could pos- 
sibly be more fundamentally tied to the stellar mass of 
the host galaxy rather than its optical color; determin- 
ing which is more fundamental is beyond the scope of 
this paper. However, the enhanced clustering amplitude 
found here for AGN compared to galaxies with the same 
color and magnitude distribution can not be due to a dif- 
ference in stellar mass between the populations; ie, the 
clustering difference can not be due to AGN being hosted 
in more massive galaxies. The stellar mass of DEEP2 
galaxies is tightly correlated with the (U — B) color an d 
Mb magnitude (see Figure 10 in lWeiner etall ((200l '). 
such that the matched galaxy samples have similar stellar 
mass distributions as the AGN. 

We further compare the clustering of the AGN in red 
and blue galaxies to the full red and blue galaxy sam- 
ples and also to red and blue galaxy samples with the 
same color and magnitude distribution ('matched' red 
and blue samples) in Figure [121 The relative biases for 
these samples are given in Table 2. There are only minor 
differences between the biases for matched galaxy sam- 
ples and the full red or blue galaxy samples. We find 
that AGN in red host galaxies are more clustered than 
both all red galaxies and matched red galaxies at the 
2.2(7 level when combining measurements from all scales, 
and on large scales they are not significantly more clus- 
tered. The relative bias is similar to that from comparing 
all AGN to all matched galaxies, but the error bars are 
larger on the red subsamples and the significance is lower. 
Comparing the clustering of AGN in blue host galaxies 
with matched blue galaxies, we do not find a significant 
difference in their clustering properties, though again the 
relative bias is consistent with that for all AGN and all 
matched galaxies; the error bars for the blue subsamples 
are simply too large to say if there is a significant dif- 
ference in clustering or not. The difference that we find 
between the clustering of the main AGN sample and all 
matched galaxies may therefore either be dominated by 
a difference in the clustering of red galaxies with and 
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without AGN or be due to both red and blue galaxies. 
6. DISCUSSION 

We find here that non-quasar X-ray AGN, taken as 
an ensemble, have a similar clustering amplitude as red, 
quiescent galaxies at z ~ 1. Our resul t s con firm and ex- 
tend the findings of lGeorgakakis et "aTl ()2007D . who found 
that a smaller sample of X-ray AGN in AEGIS had envi- 
ronments similar to those of red galaxies. Not all of the 
X-ray AGN are hosted by red galaxies, however; roughly 
half are in galaxies whose optical color is blu e. Visual in- 
spect ion of ACS imaging in the AEGIS field (jPavis et all 
l2007t ) confirms that the optical light in the X-ray AGN 
host galaxies is not dominated by the AGN itsel f, but b y 
galaxy light (see also Figure 2 of INandra et"aI1 (12001 ): 
the exception to this is the 14% of X-ray AGN that we 
define as quasars, where the optical light is from the AGN 
itself and is very bright and blue. We find that the X-ray 
AGN in red host galaxies are more clustered than those 
in blue host galaxies, and that when compared to galaxy 
samples with the same color and luminosity distribution 
(and hence, stellar mass), those galaxies with X-ray AGN 
are somewhat more clustered. 

6.1. Comparison to Other X-ray Clustering Studies 

Figure [13] compares the correlation scale length found 
for various AGN and galaxy samples at z ~ — 1. We 
find here that X-ray AGN at a median z = 0.94 with 
a median luminosity of log — 42.8 erg s^^ have a 
correlation length of tq = 5.95 ± 0.90 Mpc and 
7 = 1.66 ± 0.22. This clustering s c a,le length is simila r 
to t hat found by iGilli et all (|2005D . lYang et all (|2006D . 
and lGilh et cdl ('2009') for the CDFN and CL ASXS and 
COSM OS fields and differs from that found bv lGilh et all 
()2005[ ) for the CDFS field; these are the only other stud- 
ies at z ~ 1 that use spectroscopic redshifts to mea- 
sure the correlation f unctio n of X-ray AGN. At lower 
redshift. lMullis et al.l ()2004f ) find a consistent correlation 
length of To = 7.4 ± 1.9 Mpc for a fixed slope of 
7 = 1.8; however, a comparison with this result is com- 
plicated due to the brighter nature of the sources (me- 
dian log ix = 44 erg s~^) and larger measurement scales 
(rp = 5 — 100 Mpc). The correlation slope that 
we find, 7 = 1.66 dz 0.22, is cons i stent with that found 
bvlGilh et all 0i()^ . \Yai^set^ ^^^^. and lGilh et all 
(|2009t ). which are measured on similar scales to those 
used here. 

Our r esults are in mod est disagreement with 

those of iMontero-Dorta et al.l (|2008D . who study the 
environment-dependence within the red sequence of var- 
ious AGN classes in DEEP2, including optically-selected 
LINERs, Seyferts, and X-ray AGN, and find that while 
LINERs in red host galaxies are more likely to be in over- 
dense regions, Seyferts and X-ray AGN do not have any 
environment-dependence relative to red galaxies with the 
same magnitude distribution as the AGN population. In 
contrast, we find here that X-ray AGN in red galaxies are 
biased relative to a galaxy population matched in color 
and magnitude at the 2a level. However, only 36 X-ray 
sources were used in that study, which may account for 
the discrepancy with our findings here. 

Measuring the cross-correlation function of X-ray 
sources with galaxies, as is done here, has several bene- 
fits. One is that the spatial selection function and red- 



shift completeness for one of the samples, either the AGN 
or the galaxies, does not need to be known, as long as the 
selection function of the other sample is well character- 
ized. By using the cross-correlation of the X-ray sources 
with galaxies, we do not need to include the selection 
function for the X-ray AGN. When measuring the auto- 
correlation function of X-ray sources, one must carefully 
model both the spatial selection function, including the 
varying X-ray sensitivity across the field, and the spec- 
troscopic completeness, so as not to bias the results. All 
auto-correlation function measurements from Chandra, 
XMM, and ROSAT are potentially affected by the sensi- 
tivity varying across the field and must carefully address 
this systematic effect. 

The other benefit of using the cross-correlation func- 
tion with galaxies in the same field is one can sensi- 
tively probe the relative bias between AGN and galax- 
ies as a function of their properties, and to first order 
cosmic variance will cancel in the relative bias. Cross- 
correlation functions provide accurate relative clustering 
measuremen ts and do not in troduce systematic effects; 
as shown in ICoil et al.l (|2008l ). the relative bias inferred 
from cross-correlation measurements is consistent with 
that derived fro m the ratio o f the a uto-correlation func- 
tions. Recentlv. iHickox et al.l (|2009f ) measured the cross- 
correlation function of X-ray AGN with galaxies in the 
AGES survey (ref) at z = 0.5 to infer an absolute bias 
of6 = 1.34 ±0.16 and a mean dark matter halo mass of 

1 X IQ-^^h^^ Mq, similar to the results found here. 

Several studies have investigated the dependence of X- 
ray AGN clustering on the hardness or hardness ratio of 
the source. Within the errors, we find no significant dif- 
ference bet ween AGN s ample s div ided accord i ng to hard- 
ness ratio. iGilli et al.l (|2005( l and lGilh etall (|2009| ) also 
find no significant difference in the clustering strength of 
soft and hard X- ray sources i n the CDFN, CDFS, and 
COSMOS fields. lYang et al.l (|2003[) detect a higher an- 
gular clustering signal for hard sources compared to soft 
sources, though without a well-known redshift distribu- 
tion for the hard and soft samples this can not be inter- 
preted as a significant difference in the intrinsic cluster- 
ing properties. If, for example, the harder sources had 
a narrower redshift distribution, that would account for 
the different angular clustering signals. 

We also find no significant dependence of clustering 
amplit ude on X-ray luminosity within the range probed 
here. iPlionis et al.l (|200 8) claim to find a trend with 
brighter X-ray sources having higher clustering ampli- 
tudes, but these results rely on angular clustering mea- 
sures without known redshift distributions, and the au- 
thors infer extremely large and impl ausible correlation 
lengths at z ~ 1 of ~ 10-20 h^^ Mpc. lYang et all (|2006f ) 
detect a 'weak' correlation between clustering amplitude 
and X-ray luminosity when comparing thei r X-ray AGN 
clust ering results with that of 2dF quasars (jCroom et al.l 
|2004| ). where X-ray luminosities for the quasars are esti- 
mated from their bolometric luminosities. They find that 
quasars are somewhat more clustered than X-ray AGN, 
though at < 2a significance. This differs with what we 
find below (and the results shown in Figure [T3)) . when 
comparing the clustering of X-ray AGN and quasars in 
the sa me volume at t he sa me redshift. The z = 0.2 re- 
suhs of iMullis etUn ()2004f) for brighter X-ray AGN do 
indicate a slightly higher correlation length than those at 
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Fig. 13. — Clustering results for various AGN and galaxy sur- 
veys with < 2 < 1. The comoving correlation scale length, rp, 
is plotted for different samples as a function of redshift. Shown in 
black with dashed line error bars are X-ray AGN clustering results 
from MuUis ct al. (2004) (open bo x). IGilli et a.1.1 ('2005') (triangle), 
Euig ct al. (2006) (plus sign), and lGi h et al.l l| 2009) (X sign). Re- 
sults for X-ray AGN from this work are shown with a green filled 
circle. Results for red and blue galaxi es are shown as red and blue 
markers for SDSS galaxies from Zehavi e t al.l II2005I) (plus sign) 
and DEEP2 galaxies from Coil ct al. (200q) (open diamond). At 
z = 0.1 galaxy clustering results from 2dF are similar to those 
shown for SDSS. Orange mark ers indicate measures of optically- 
bright quasar clustering from Porciani. Maglio cchetti. fc Norberal 
II2004I ) (orange open box) and Coil et al. (2007) (orange X sign). 
The clustering of X-ray AGN found in this paper show that non- 
quasar X-ray AGN cluster similarly to red galaxies at z ~ 1 and 
are more clustered than blue galaxies. 

z ~ 1, but that may be due to evolution and growth of 
structure, not intrinsic luminosity-dependent clustering. 

The lack of a convincing relation between X-ray lumi- 
nosity and clustering strength is important for theoretical 
models, as the X-ray luminosity is related to the black 
hole accretion rate for a fixed black hole mass. This may 
indicate then that there is not a strong correlation be- 
tween accretion rate and host halo mass, i.e., between 
how much gas is accreted onto a black hole at a given 
time and how much gas potentially exists in the reser- 
voir of the halo. 

6.2. Comparison to Optical AGN Clustering Studies 

A comparison of the clustering properties of different 
kinds of AGN can lead to a clearer understanding of their 
respective host galaxies, fueling mechanisms and AGN 
lifetimes. The clustering of X-ray AGN was not previ- 
ously known to high enough precision to make a detailed 
comparison with other AGN samples. Among optically- 
selected AGN, clustering studies have been performed 
both for quasars at intermediate and high redshift and 
lower-luminosity Seyferts at low redshift; 

Using the same techniques as here, ICoil et aD ()2007l ) 
find that spectroscopic broad-line SDSS quasars in the 
DEEP 2 volume at0.7<z<1.4 cluster like blue galax- 
ies rather than red galaxies (at 2a significance), with 
ro = 3.1 ± 0.6 h^^ Mpc. Comparing the scale length 
found here for no n-quasar X-ray AGN with that found by 
ICoil et al.l (120071 ) for quasars, the X-ray AGN are more 
clustered than quasars at the 2.6a level. This comparison 
should be quite robust, as the same sample of galaxies 
from the DEEP2 survey is used as a tracer of large-scale 



structure in both analyses, and therefore the survey se- 
lection, and design are identical, as are the scales probed. 
This comparison suggests that high-accretion AGN, ob- 
served as optically bright quasars, are likely hosted in 
star-forming galaxies, while X-ray AGN are hosted both 
by massive bl ue galaxies and red ga laxies, as seen in Fig- 
ure [H] (see also iNandra et al.l (|2007f) ). reflecting a general 
correlation between the amount of star formation in a 
galaxy a nd the accretion r a te of its ce ntral black hole 
see e .g., iKauffmann et al.l (200'f) and I Silverman et al.l 
20081 ) for other observational signatures of this general 
correlation). As both quasar fueling and star formation 
require the presence of cold gas, it is reasonable that they 
should be found in similar environments. // optically- 
bright, high-accretion quasars and lower-accretion X-ray 
AGN are similar objects at different evolutionary stages, 
then these results are consistent with a scenario in which 
a quasar resides in a star forming galaxy in the blue cloud 
before having its star formation quenched and moving to 
the red sequence with a lower luminosity X-ray AGN. 
The timescale between the quasar and lower-accretion 
X-ray AGN stages would have to be on the order of sev- 
eral Gyr, however, to allow enough time for the clustering 
amplitude to increase due to gravity. 

Other studies of q uasars at z ^ 1 (e.g., 
iPorciani. Magliocchetti. fc Norberj 12004 iMvers et all 
120061: Ida Angela et all 120051) find somewhat higher 
corr elation lengths ( consistent at the 2a level with 
the ICoil et al.l (|2007l ) results) using the quasar auto- 
correlation function, though usually with larger er- 
ror bars due to the low number density of quasars. 
The most constraining of thes e resu lts is from 
IPorciani. Magliocchetti. fc Norber j (|2004j ). using the 
2QZ survey, who find ro = 4.7±0.7 h'^ Mpc at z = 1.06 
(shown as an ora nge open box in F igurell3 [) . This result 
is consistent with ICoil et all (j2007( ) but allows a correla- 
tion length that matches either red or blue galaxies, given 
the error bars. This is generally true of other quasar clus- 
tering results at this redshift. Quasar clustering studies 
generally find a correlation slope, 7, of ~ 1.7, very sim- 
ilar to that of blu e galaxies but lo wer than that of red 
galaxies at z ~ 1 (|Coil et al.ll2008( ). which implies that 
quasars can not predominantly be hosted by red galaxies. 

While a comparison of the quasar and AGN cluster- 
ing results from DEEP2 point to AGN possibly being 
more clustered than quasars, a more precise measure- 
ment of the quasar clustering amplitude is needed to 
more firmly identify the host galaxy types and dark mat- 
ter halo masses of quasars. The fundamental limitation 
in our ability to more precisely measure the quasar corre- 
lation function is the low number density of quasars; this 
makes cross-correlations of quasar samples with more 
abundant galaxy samples particularly effective. Future 
precise measurements of the quasar correlation function 
will likely come from wide-area galaxy redshift surveys 
at intermediate and high redshift. 

6.3. Black Hole Masses for Quasars and X-ray AGN 

The 2a difference in clustering amplitude measured at 
z ~ 1 for quasars and X-ray AGN, if real, would reveal 
a difference in the mean host dark matter halo mass of 
these two populations, which could reflect a difference in 
the mean black hole masses for quasars and X-ray AGN. 
To investigate this we estimate black hole masses for the 
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quasars used in ICoil et all ()2007[ ) and the X-ray AGN 
used in this study. 

To estimate black hole masses for t he br oad-line SDSS 
quasar population used in lCoil et al.l ()2007f ). we use prop- 
erties measured in the SDSS spectra, in particular the Mg 
II 2800 A linewidth and UV flux, which are available for 
all of the qu asars used in that s t udy. The following es- 
timator from lMcLure k, Dunlod (|2004f ) is used, which is 
calibrated at low redshift using reverberation mapping: 
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where L3000 is the restframe continuum lu minosity at 
3000 A . We find that the SDSS quasars in the lCoil et aH 
(|2007t l sample have a median black hole mass of Mbh = 
3.0 X IO^Mq and a mean mass of Ubh = 4.3 x lO^M©, 
with a Icr standard deviation of 3.4 x 10^ Mq, which 
agrees w ell with the masses of SP SS quasars at z ~ 1 
found bv lMcLure fc Dunlopl (f200l ). There may be sys- 
tematic biases associated with this method, particularly 
given that we are applying it at a different redshift than 
it was calibrated, but the estimated masses should be 
within a factor of ~2-3 of the true masses. 

Black hole mass estimates for the X-ray AGN sam- 
ple used in t his pa per are presented and discussed in 
iBundy et al.l (|2008l ). where the following relation be- 
tween black_hole_mass and bulge infrared luminosity is 
used (|Grahan]||2007l) : 



log 



M 



BH 



Mr. 



-0.37(Mk, -I- 24) + 8.29. (9) 



K-band photometry for the Chandra sources is 
from a Palomar/W IRC survey of the DEEP 2 fields 
(jBundv et al.ll200^ . To estimate the black hole masses 
we assume no evolution in the black hole-bulge mass rela- 
tion with redshift. We also assume that the host galaxies 
of the X-ray AGN are bulge-dominated and that the bulk 
of the K-band light is from the bulge. This is a reasonable 
assumption, as HST imaging shows that m ost of the X- 
ray A GN are in bulge-dominated galaxies (jPierce et al.l 
|2007[ ). There is a systematic effect, however, that would 
drive the black hole mass estimates lower as less of the 
K-band light is in the bulge. Here again we expect these 
black hole masses to be within a factor of ~2-3 of the 
true masses. We find that the median black hole mass 
for the X-ray AGN is Mbh = 4.7 x IO^Mq and the mean 
is Mbh = 6.2 x 10^ Af©, with a Icr standard deviation of 
5.8 X 10** Mg. 

The inferred X-ray AGN black hole masses are compa- 
rable to the quasar black hole masses, which could imply 
that they are in similar stellar mass galaxies. It may be 
surprising then if the non-quasar X-ray AGN are more 
clustered than the quasars. If true, this could indicate 
that quasars are more prevalent in lower mass halos be- 
cause those halos have more fuel available to power the 
b lack hole. 

iFerrares^ ()2002f ) finds a scaling between black hole 
mass and host dark matter halo mass at low redshift: 

Mbh/10«Mq = 0.67 (Mhaio/10i2MQ)i-82^ (10) 



where results from ISeliakI (|2002f ) have been used for 
halo profile fits derived from weak lensing measurements. 



Given our rough black hole mass estimates above, this 
relation implies that both the SDSS quasars and non- 
quasar X-ray AGN should be in ~ 3 x lO^^M© ha- 
los, if there is no evolution in t his relation to z ^ 1. 
As presented in ICoil et al.l |2008) , from clustering mea- 
surements the minimum dark matter halo mass for blue 
galaxies at z 1 is Mmin ^ 5 x lO^^h^^ Mq and for 
red galaxies is Mmin ^ 2 x 10^^ /i~^Mq. The corre- 
sponding mean halo masses for blue and red galaxies are 
M 1.8 X IO^^Mq and M ~ 5 x lO^^M©, respectively, 
for h — 0.7. Assuming that the scaling law between black 
hole mass and dark matter halo host mass does not evolve 
to z ~ 1, this would imply that quasars and non-quasar 
X-ray AGN reside either primarily in blue galaxies or in 
a mix of red and blue galaxies, which is consistent with 
the location of the X-ray AGN in the color-magnitude 
diagram. 

6.4. X-ray AGN and Green Galaxies 

The clustering of 'green' galaxies in the minimum of 
the optical color bimodality at z ~ 1 has been measured 
by Coil ct al. ( 2008). A comparison of the green galaxies 
- which may be in transition from the blue, star forming 
galaxy population to the red, quiescent galaxy popula- 
tion - to our measurements for the X-ray AGN shows 
that both populations have a high clustering amplitude 
on large scales and therefore similar values of the correla- 
tion length (ro = 5.17 ± 0.42 Mpc for green galaxies 
with Mb < —20), comparable to red galaxies. Unlike red 
galaxies, however, green galaxies have a relatively shal- 
low clustering slope, 7 = 1.59 ± 0.08, similar to that for 
blue, star-forming galaxies. The X-ray AGN clustering 
slope is not well constrained here and is consistent with 
the slope of either red or blue galaxies. 

The location of the X-ray AGN in the optical color- 
magnitude diagram (see Figure [5]) indicates that there is 
likely to be overlap between the X-ray AGN and green 
galaxy populations. Their clustering properties reinforce 
this; both populations reside in massive halos. The shal- 
low slope of the green galaxy correlation function likely 
reflects a greater tendency for green galaxies to be found 
near the outskirts of dark matter halos, rather than their 
centers, compared to red galaxies. The same may be true 
for the X-ray AGN, though the apparent rise of Wp{rp) 
on small scales might indicate that X-ray AGN are ac- 
tually more likely to be central galaxies. However, given 
the large uncertainties in Wp{rp) in the innermost bins 
we can not distinguish between these scenarios here. 

6.5. X-ray AGN in Groups 

The higher clustering amplitude, found here at a sig- 
nificance of 2.8 CT, for X-ray AGN compared to galaxies 
of the same color and magnitude distribution suggests 
that a higher fraction of X-ray AGN host galaxies are in 
groups than all galaxies of the same color and magnitude, 
as galaxies in groups have a higher clustering amplitude 
(|Coil et al.l2006bHYang et al.ll2005D . As discussed in Sec- 
tion 5.3, it can not be due to a difference in the stellar 
mass distribution of the AGN host galaxies compared to 
the matched galaxy sample, as they have similar stellar 
masses (as inferred from their host galaxy colors and lu- 
minosities). For a given massive galaxy sample, whether 
star forming or quiescent, it is the galaxies that reside in 



16 



more massive halos, and therefore are more likely to be 
i n groups, that host X-ray A GN. 

iGeorgakakis et all ()2008b[ ) measure the fraction of X- 
ray AGN at z ~ 1 in groups in AEGIS using the same 
Chandra and DEEP2 gala xy data used her e. Using 
a DEEP2 group catalog (Gerke et all l2005f l that has 
been tested and calibrated with simulations and mock 
galaxy catalogs, they find that when compared to the 
full DEEP2 galaxy population X-ray AGN are preferen- 
tially found in groups, but when compared to galaxies of 
the same color and luminosity distribution of the AGN 
hosts this difference drops to the < 2a level. This is al- 
most certainly consistent with what we find here, as the 
higher clustering amplitude found for X-ray AGN host 
galaxies when compared to all galaxies of the same color 
and magnitude distribution does not imply that X -ray 
AGN are in groups identified bv lGerke et al.l (|2005l ). but 
rather that they generally reside in higher mass host ha- 
los which have a stronger likelihood of being in a group. 
Additiona lly, random erro rs are larger in the study of 
[Georgakakis et al.l (|2008bf) . due to small number statis- 
tics for the group catalog, such that what is significantly 
detected here may be seen at lower significance in the 
group study. 

6.6. The Fueling of X-ray AGN 

We find that non-QSO X-ray AGN host galaxies reside 
in more massive dark matter halos than the full popula- 
tion of galaxies of the same color and magnitude distribu- 
tion. This likely reflects a higher incidence of X-ray AGN 
in galaxy groups than in isolated galaxies. Residing in a 
galaxy group may lead to accretion onto a central black 
hole, if galaxy mergers (which are more likely to occur in 
groups) genera lly supply gas to the central AGN (e.g. , 
Barnes & Hernquist 1992; 'Mihos fc Hcrnquistl Il996t 
SEringel, Di Mattco, & Hernquist 2005b). Alternatively, 
as X-ray AGN favor host dark matter halos that are on 
the massive end of the distribution for given observed 
galaxy properties, their presence in these halos may not 
reflect interactions with the environment, but rather a 
tendency for more massive black holes in otherwise sim- 
ilar galaxies to be found in more massive halos, possibly 
due to the presence of more gas that could fuel the black 
hole. For X-ray AGN in red host galaxies, there may be 
hot gas on large scales which could cool on small scales 
near the black hole to provide fuel. 

The host galaxy morphologies o f ^^60 X-ray AGN in 
AEGIS is studied bv lPierce et all (pOOl ( see also Geor- 
gakakis et al. in preparation), who find that roughly half 
of the host galaxies are elliptical or bulge-dominated, and 
that while the fraction that are classified as mergers is 
higher than for the field (at the 2a level), this accounts 
for only ^20% of the population. It is therefore not clear 
that X-ray AGN are necessarily fueled by on-going or re- 
cent major mergers; they may be fueled by minor mergers 
or, in the case of those in red host galaxies, they may be 
fueled by small amounts of hot gas that remain long after 
a major merger. Indeed, roughly half of the X-ray AGN 
lie in red host galaxies that may have been on the red se- 
quence for a long time and possibly not have undergone a 
major merger for many Gyr. To reconcile our results with 
a scenario in which X-ray AGN are fueled primarily by 
major mergers, one would require that those in red host 
galaxies had recently migrated to the red sequence (e.g., 



Saliin et al.l 120071; ISchawinski et all 120071; IShen et all 
2007HG eorgakakis et al. I l2008at iKocevski et al. II2Q08D. 

As discussed above, the results presen ted here, when 
combined with those of ICoil et al.l (j2007f) . are consistent 
with an evolutionary sequence in which bright quasars 
are ignited while being hosted by star forming galaxies 
and later evolve to be lower-luminosity AGN, detected 
in X-rays but not as broad-line optical sources, hosted 
by massive galaxies in either the bright or red end of 
the blue cloud or on the red sequence. This evolution- 
ary sequence is not a unique explanation of the observa- 
tions, i.e.. X-ray AGN may be a different population than 
quasar relics, but if some fraction of blue star-forming 
galaxies have their star formation quenched and migrate 
to the red sequence, then it is plausible that some quasars 
may eventually become lower-luminosity AGN detected 
not by their optical light but by their X-ray emission. In 
this picture the AGN type could be due to gas availabil- 
ity in the host galaxy, which may be related to host hal o 
mass (e.g., ICroton et all [20061 ; |DekeL& Birnboim"2006'). 

The results presented here and in Coil ct al. (2007|) 
ar e generally consiste nt with the cartoon model shown 
in iHickox et al.l (|2009( l for how galaxies and AGN evolve 
with time. In this cartoon model optically-bright quasars 
are hosted by on-going disk galaxy mergers and immedi- 
ately precede an optically-faint X-ray AGN phase, which 
evolves into an early-type galaxy. 

6.7. Comparison to Theoretical Models of ACN 
Evolution 

Theoretical models of AGN formation and evolution 
can yield measurably different predictions for the lo- 
cal environments and cluste ring properties of AGN. 
iKauffmann fc Haehneltl (|2002l ) use a semi-analytic model 
in which AGN are fueled by galaxy mergers, where the 
peak AGN luminosity depends on the mass of gas ac- 
creted by the black hole, which in turn depends on the 
host halo mass. This leads to a natural prediction that 
brighter AGN reside in more massive halos, such that the 
AGN clustering amplitude should be strongly luminosity- 
dependent. This model is not supported by observa- 
tions, which generally show a lack of a strong correla- 
tion between AGN clustering amplit ude and luminosity , 
except at the very bright end (e.g., ICroom et al.l 120021; 
Porciani fc Norbergl [20061 ; iMvers etall 120071 ; IShen et alJ 
20081). 

Ho pkins et al.l (|2005f ) present an alternative model 
in which bright and faint AGN are in similar physi- 
cal systems but are in different stages of their life cy- 
cles. This model predicts that faint and bright AGN 
should reside in similar-mass dark matter halos and that 
quasar clustering shoul d depend only weakly on lumi- 
nosity (|Lidz et al.ll2006f ). That general prediction agrees 
well qualitatively with observations of quasar clustering, 
but the model also predicts that lower luminosity AGN 
should have an equal or lower clustering amplitude than 
bright AGN, which is not well-supported when compar- 
ing our results here fo r non-quas a r X-r ay AGN with the 
results for quasars in ICoil et all (I2007D. In the specific 
model presented bv lHopkins et alF i 20081 ). AGN are de- 
tected in X-rays while obscured by dust just after a ma- 
jor merger, before undergoing an optically-bright quasar 
phase a short time later; this picture is also not well- 
supported by our results. 
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The semi-analytic model of ICroton et al.l (|2006f ) com- 
bines a pr escription of merger-driven bla ck hole growth 
similar to iKauffmann k, Haehneltl (|200Cl| ) with an inde- 
pendent mode for hot gas accretion in large halos which 
accounts for the fueling of lower-luminosity AGN in mas- 
sive halos. This model assumes that some fraction of cold 
gas must be present to trigger a bright quasar phase dur- 
ing a galaxy merger; the black hole itself must also be 
massive such that the luminosity remains sub-Eddington. 
Quasars, in this model, can be strongly clustered at high 
redshift when cold gas fractions are presumably high, but 
as their host dark matter halos grow cooling becomes 
more inefficient as the virial temperature of the halo in- 
creases, and the cold gas supply is suppressed above a 
given threshold mass. Quasars will only be found in ha- 
los below the threshold mass, and thus relatively gas- 
free red galaxies are not expected to shine as quasars; 
only those mergers that occur in lower mass halos, pre- 
sumably outside of group environments at 2: < 1, will 
contain sufficient cold gas to fuel a quasar. This model 
therefore predicts that quasars should cluster similarly to 
massive star-forming galaxies at a given redshift, where 
both the star formation and quasar activity are fueled 
by cold gas, while massive halos host lower-luminosity 
AGN that are fueled (possibly relatively inefficiently) by 
hot gas (which may cool on small scales near the black 
hole). X-ray AGN, in this picture, therefore should clus- 
ter like quasars did at an earlier epoch. There are no 
quantitative predictions from this model for the cluster- 
ing of, say, the full X-ray AGN population versus bright 
quasars, but the model is qualitatively consistent with 
our results and the relati vely higher cluster i ng seen for 
quasars at high redshift (|Myers et al.l l2007t IShen et all 
iCroom et a Il2n05l) . 

Thacker et al. r(l2008D present a model for quasar fu- 



eling and feedback that quantitatively matches observa- 
tions of quasar and X-ray clustering reasonably well. In 
their model bright AGN (both quasars and bright X-ray 
AGN) are the result of mergers and are augmented by 
feedback from AGN outflows. They do not present pre- 
dictions for z < 1.2, but their model matches well the 
observed clustering of quasars at z — 1.5 — 2. Their 
prediction for the correlation length of X-ray AGN at 
z — 1.75 with a luminosity Lx = 3.2 x 10^^ ergs"^ is 3.6 
Mpc, significantly lower than what we find here at 
z — 0.9. However, as it is not clear that the bulk of the 
X-ray AGN studied here are the result of major mergers, 
this model may not be entirely applicable to the results 
presented here. 

7. CONCLUSIONS 

Using the cross-correlation of galaxies and Chandra- 
detected AGN in the AEGIS survey, we are able to de- 
termine the clustering of X-ray AGN at z 1 much 
more precisely than studies that measure the cluster- 
ing of AGN using the auto-correlation function of AGN 
alone. We find that non-quasar X-ray AGN with a me- 
dian log Lx ~ 43 erg have a clustering scale length 
ofro = 5.95 ±0.90/1-1 Mpc and slope of 7= 1.66 ±0.22, 
with a bias that corresponds to a minimum dark matter 
halo mass of M™„ = 5 (+5/-3) x IQI^/i-IMq. The clus- 
tering amplitude of X-ray AGN is consistent with that 
of red and 'green' galaxies at the same redshift and is 



greater than that of blue galaxies, even though roughly 
half of the X-ray AGN are hosted by blue, star-forming 
galaxies. We do not find a significant difference in clus- 
tering amplitude within the AGN population when di- 
vided into subsamples based on on optical luminosity. 
X-ray luminosity, or hardness ratio, within the ranges 
probed here. 

The clustering strength of X-ray AGN is primarily de- 
termined by the host galaxy color; there is a smaller but 
significant effect due to a difference in clustering am- 
plitude for galaxies that host X-ray AGN compared to 
matched samples of galaxies with the same color and lu- 
minosity distribution. This implies that galaxies that 
host X-ray AGN arc in more massive dark matter halos 
than those that do not, for fixed host galaxy properties; 
i.e., galaxies that host X-ray AGN are more likely to be 
in groups when compared to galaxies of the same color 
and luminosity distribution. 

By comparing to previous results on the clustering of 
optically-selected quasars in the DEEP2 fields, we find 
that non-quasar X-ray AGN are more clustered tlian 
optically-selected quasars at the same redshift at the 2.6cr 
level. Our results are consistent with quasars occurring 
in blue star-forming galaxies which later settle onto the 
red sequence and potentially host lower-luminosity X-ray 
AGN, if they are similar objects at different evolutionary 
stages. While clustering measures of different types of 
AGN can be used to constrain their host galaxy types, 
halo masses, and fueling mechanisms, precise measure- 
ments await larger galaxy and AGN redshift surveys at 
intermediate and high redshift. 
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